MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 


THE  APPUCATION  OF  THE  CRAMER-RAO  BOUND 
TO  ESTIMATES  OF  RADAR  RETURN  TIME -OF -ARRIVAL 
FOR  SEVERAL  TARGET  CONFIGURATIONS 


C.  B.  CHA^G 
R.  W.  MILLER 

Group  32 


TECHNICAL  NOTE  1977-27 


23  MAY  1977 


Approved  for  public  release;  distribution  unlimited. 


LEXINGTON 


MASSACHUSETTS 


ii 


Abstract 


In  this  report,  certain  applications  of  the  Cramer-Rao 
bound  to  problems  in  radar  echo  time-of-arx’ival  estimation  are 
presented.  These  applications  encompass  point  target  echos, 
two-point  target  echos,  and  distributed  target  echos.  With 
respect  to  point  targets  and  two-point  targets,  a modified  Cramer- 
Rao  bound  is  derived  and  the  results  are  substantially  tighter 
than  those  of  the  conventional  Cramer-Rao  bound.  In  applying  the 
Cramer-Rao  bound  to  distributed  targets,  a sampled-data  formulation 
is  presented  and  a numerical  example  corresponding  to  a wake 
parameter  estimation  problem  is  given  to  illustrate  the  results. 
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Introduction 


The  familiar  Cramer-Rao  bound  on  the  variance  of  estimates 
is  particularly  useful  since  it  is  frequently  easy  to  calculate, 
and  is  generally  a tight  bound  for  a wide  class  of  estimators 
when  the  signal-to-noise  ratio  is  high  (References  ll]-(4]). 

Although  there  are  c*-her  bounds  also  tight  for  small  signal-to-noise 
ratios,  they  are  usually  more  difficult  to  obtain  (References  II] 
-{6]).  For  this  reason,  the  Cramer-Rao  bound  is  often  used  in  radar 
applications . 

In  estimating  unknown  (deterministic  or  random)  parameters 
from  noisy  measurements,  there  often  exist  nuisance  parameters 
(deterministic  or  random) , In  radar  applications,  the  parameter 
of  interest  may  be  the  radar  return  time-of-arrival,  or  target 
range,  while  the  nuisance  parameters  may  be  signal  amplitude  and/or 
phase.  Although  we  often  do  not  care  to  estimate  these  nuisance 
parameters,  the  estimates  of  interesting  pareuneters  are  nevertheless 
dependent  on  them.  In  the  usual  application  of  the  Cramer-Rao  bound, 
the  marginal  density  of  the  measurement  conditioned  on  the  interest- 
ing parameter  is  used  so  that  the  dependence  on  the  nuisance  para- 
meter is  never  explicitly  shown.  This  may  sometimes  result  in 
optimistic  or  even  trivial  bounds.  In  this  report,  we  obtain  a 
modified  bound  by  treating  the  nuisance  parameters  differently. 

These  modified  bounds  apply  to  a different,  more  restricted  class 
of  estimators  than  the  conventional  Cramer-Rao  bound,  but  are 
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useful  in  several  radar  applications  where  bounds  tighter  than  the 
ordinary  Cramer-Rao  bound  are  obtained. 

The  application  of  the  Cra*ner-Rao  bound  (ordinary  and  modified) 
for  estimates  of  parameters  of  a random  signal  is,  however,  usually 
difficult  to  obtain.  Random  signals  often  occur  in  practice,  for 
example,  the  radar  return  from  the  wake  trailing  a high  speed  re- 
entry vehicle.  Returns  of  this  kind  are  often  thought  of  consisting 
of  a large  number  of  radar  scatterers  with  their  relative  amplitude, 
phase,  and  location  being  randomly  distributed.  They  are,  therefore, 
sometimes  referred  to  as  distributed  targets.  It  was  shown  (Reference 
[ 2 1) that  the  Craraer-Rao  bound  for  a random  signal,  or  a distributed 
target,  may  be  obtained  via  the  solution  of  a rather  complicated 
integral  equation.  Except  for  a few  special  cases,  the  solution 
can  only  be  obtained  by  using  approximations.  In  this  note,  we 
consider  the  calculation  of  the  Cramer-Rao  bound  for  a certain 
class  of  random  signals  via  a sampled  data  approach.  The  random 
signal  considered  is  a conditional  Gaussian  process  (Reference  [2]) 
with  known  time  varying  mean  and  covariance  functions.  Since 
in  many  radar  applications,  the  data  are  often  recorded  and  pro- 
cessed in  sampled  form,  the  approach  presented  in  this  report  is 
believed  to  be  practical  and  useful.  A nxamerical  example  correspond- 
ing to  a wake  parameter  estimation  problem  is  given. 
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This  report  is  organized  as  follows.  The  modified  Cramer-Rao 
bound  is  described  in  Section  2.  Its  application  to  several  point 
targets  is  given  in  Section  3.  The  application  of  the  modified 
bound  to  the  estimate  of  time  separation  of  two  closely  spaced 
signals  is  presented  in  Section  4.  In  Section  5 we  present  the 
sampled  data  approach  to  the  calculation  of  the  Cramer-Rao  bound 
for  distributed  targets.  A numerical  example  is  included  to 
illustrate  the  result.  Two  appendices  are  attached.  In  Appendix  A 
the  laborious  derivation  required  to  obtain  the  final  results  of 
Section  4 is  given.  There  are  at  least  two  limiting  cases  for  dis- 
tributed targets.  Under  proper  conditions,  a distributed  target 
converges  to  a known  target  or  a Rayleigh  target.  Derivations  to 
explore  these  relations  are  given  in  the  Appendix  B. 

2.  A Modified  Cramer-Rao  Bound 

In  this  section,  application  of  the  Cramer-Rao  bound  to  the 
case  wherein  there  are  nuisance  parameters  is  discussed.  By  treating 
nuisance  pareuneters  in  an  unusual  way,  we  develop  a modified  bound. 

We  shall,  for  simplicity,  discuss  only  the  scalar  estimation 
problem  although  results  are  easily  extended  to  the  multiple 
parameter  estimation  problem.  Two  cases  are  considered. 

CASE  1.  Nonrandom  Parameter 

Let  a be  an  unknown  parcuneter  which  is  to  be  estimated  based 
on  a set  of  observations  r.  Let  the  joint  probability  density 
function  (PDF)  be  known,  and  equal  to 


(2.1) 


(R/a) 


The  Cramer-Rao  bound  is  a bound  on  the  variance  of  any  unbiased 
estimator  of  a:  a(r).  It  is: 
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a 
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(R/a) 
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(R/a) 
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=C(a) 


(2.2) 


If  the  observation  depends  not  only  on  *’d'  but  also 
"nuisance"  parameter,  0,  then  the  bound  also  depends  on 
of  6.  If  0 is  random,  with  a known  PDF,  then  the  bound 
using 


on  a 

the  value 
is  calculated 


p^(R/a) 


Pj.(R/0,a)  Pg(0/a)  de 


(2.3) 


The  bound  then  applies  to  estimators  which  are  unbiased  in  a global 
sense,  but  not  necessarily  for  each  value  of  0.  This  leads  to  the 
following  alternative  procedure: 


Let  a(r)  be  an  estimator  of  the  poirameter  a which  is  unbiased 
for  every  value  of  the  nuisance  parameter  0.  Then 


> Eg  {C(a/9) } 


(2.4) 
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where  C(a/0)  is  the  Cramer-Rao  bound  calculated  for  each  value  of 


6,  and  Eq  denotes  expectation  over  the  random  variable  0. 

y 


The  above  result  is  obvious.  For  any  estimator  which  is 


>2,, 


unbiased  for  all  0,  its  variance  a"  10)  for  each  particular  0 ii 

cl 


bounded  by: 


.^2 


a:  (0)  > C(a/0) 

a * 


where  Cl  a/0)  is  calculated  using 
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C(a/0) 


-1 


4 
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lnp^(R/0,a) 


11“ 


Since  the  inequality  holds  for  all  0 and  since  the  estimate  is 
unbiased  for  all  6,  it  follows  that: 


a"'  = E, 


{Ve[ 


a(r)-a 


II 


= (0)i  > E^  / C(a/0)  \ 


CASE  2. 


Random  Parameter 


Let  a be  a random  parameter  which  is  to  be  estimated  based 
on  a set  of  observations  r.  Let  the  joint  PDF  of  r and  a be  known, 
and  denoted  by 


P_  ^ (R,A)  (2.8) 

•L  |i  d 

The  Cramer-Rao  bound  on  the  variance  of  any  unbiased  estimator  of 
a is: 


If  the  observation  depends  also  on  a parameter  0,  then  the 
bound  also  depends  on  the  value  of  0.  If  0 is  random,  with  a knovm 
PDF,  then  one  uses 


QO 

Pr.a'S.'-''  =/  Pr.a/e'S.ve)  Pg(e)  de 


(2.10) 
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The  bound  then  applies  to  estimators  which  are  unbiased  in  a global 
sense,  but  not  necessarily  to  estimators  which  are  unbiased  for 
each  value  of  6.  This  leads  t-o  the  following  alternate  procedure: 

Let  a (£)  be  an  estimator  of  the  random  variable  a which  is 
unbiased  for  every  value  of  the  random  nuisance  parameter  Then: 


0 > E.  {C(a/e)  ) 

a — 6 


'2.11) 


where 


C(a/9) 


• ! -y 


^"Pr,a/8 


R,A/8)j  I 


-1 


a/e  (E.A/9))' 


k-1 


(2.12) 


The  proof  follows  that  for  nonrandom  parameters. 


Discussion: 

The  difference  between  our  bounds  and  the  standard  Cramer-Rao 
bound  on  the  variance  of  estimators  with  a random  nuisance  parcuneter 
is  that  our  bounds  apply  to  estimators  which  are  unbiased  for  each 
value  of  the  nuisance  parameter  whereas  the  standard  bound  applies  to 
estimators  which  need  be  unbiased  only  over  the  ensemble.  For 
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instance,  in  the  nonrandoa  paromefo.'  case,  if  a(r)  denotes  the 
esti  iator,  then  the  rtanvivird  LoanJ  is  satisfied  if: 


.(a)  = f a (B’i  p (H,  'a)  dRd  - a (2.13) 

IT  ■ * ^ ~ n , 


whereas  our  boun  i is  satist  iea  iinii,-!  the  more  restricted  condition; 


^r/-  dB  = a (2.14) 

tor  all 

Certainly,  one  is  not  always  concerned  that  an  estimator  h': 
unbiased  for  all  values  ot  tne  nuisance  parameter.  Howe  sir,  there 
may  be  many  applications  where  that  property  is  desirable.  For 
instance,  in  estiraatinq  the  time-ot-arr ival  of  a radar  echo  of  random 
amplitude,  one  would  certainly  like  the  cstimato.  to  be  unbiased  for 
all  values  of  signal  amplitude.  If  the  t*..  ’^-arrival  is  a nonrandom 
parameter,  this  requirement  may  be  unreasonable  if  the  signal  amplitude 
can  be  too  small  for  reliable  detection.  However,  such  an  unbiased 
estimator  can  be  conceived  if  the  time-of-arrival  is  a random  variable, 
since  its  a priori  distribution  provides  an  unbiased  estimate  if  the 
signal  is  not  detected.  A later  sot  of  examples  treats  this  case. 

The  modified  bound  which  we  have  presented  has  a natural  inter- 
pretation. The  quantity  C(a/*)  is  the  bound  which  is  obtained  under 
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the  assumption  that  the  nuisance  parameter  0 is  known  to  the 
estimator.  The  modified  bound  is  obtained  by  simply  averaging 
this  quantity  over  0.  One  might  expect  that  a bound  obtained  under 
the  assumption  that  the  random  nuisance  parameter  is  always  known 
to  the  estimator  would  not  be  as  tight  as  one  which  does  not  make 
this  assumption.  However,  this  is  not  always  the  case;  as  our 
examples  will  demonstrate - 

There  is  another  alternative  when  a nuisance  parameter 
is  present.  That  is,  to  estimate  jointly  the  values  of  both  the 
nuisance  parameter  and  the  desired  parameter.  This  leads  to  the 
Cramer-Rao  bound  for  the  multiple  parameter  case.  To  evaluate  it 
requires  the  evaluation  of  ijnCn+l)  expectations  of  derivative  of 
logarithms  of  the  joint  PDF  followed  by  the  inversion  of  an  n by 
n Fisher  Information  matrix,  where  n is  the  number  of  parameters 
estimated.  Since,  in  practice,  the  estimator  actually  used  is 
frequently  one  which  estimates  nuisance  parameters  as  well  as  the 
desired  parameter,  e.g. , a joint  maximum  likelihood  estimator,  such 
a bound  may  be  more  representative  of  the  actual  sitUu.tion.  Un- 
fortunately, the  difficulty  in  evaluating  the  bound  increases  as 
the  number  of  nuisance  parameters  increases  and  may  become  intrac- 
table rather  quickly. 

3.  Point  Targets 

In  this  section,  we  apply  the  modified  Cramer-Rao  bound  to 
several  point  target  configurations.  Let  s(t)  denote  the  low  pass 
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transmitted  radar  signal,  the  echo  from  a point  target  is  therefore 
the  time  delayed  and  complex  modulated  version  of  s(t).  The  complex 
observation  of  the  returned  signal  is: 


r(t)  = as(t-x)  + w(t) 


where:  s(t)  is  the  complex  signal  envelope 


J |s(t)|2 


dt 


1 


(3.  1) 


(3.2) 


w(t)  is  complex  white  Gaussian  noise  with  spectral 
density  N^/2  (two  sided) 

a is  the  complex  signal  amplitude 

T is  the  unknown  Time-of-ar rival 

For  a fixed  known  value  of  a,  the  well  known  Cramer-Rao  bound 
for  estimators  of  x is  [3]: 


(E)  > C (E) 
T — 


where  E 


= energy  in  the  signal 


(3.3) 
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3 


signal  root-mean-square  bandwidth  in  radians/sec 


3 


2 


s ( t+r/2) s 


* 


(t-x/2) 


T=0 


(3.4) 


For  a LFM  radar  pulse,  the  relation  between  3 and  the  instan- 
taneous bandwidth  B,  in  Hz,  taking  into  account  the  weighting  used 
to  reduce  sidelobes,  is  3 as  1.2  B. 


In  practice,  the  complex  signal  amplitude  a may  fluctuate 
according  to  some  random  distribution.  Suppose  that  the  complex 
signal  amplitude  is  random  with  phase  uniformly  distributed  over 
[o,2tt]  and  energy  distributed  according  to  a PDF  p(E)  . Ideally, 
we  would  wish  the  time-of-arrival  estimator  to  be  unbiased  for  all 
values  of  E.  Following  the  discussion  on  the  modified  Cramer-Rao 
bound,  a bound  on  the  variance  of  such  an  estimator  is: 


a 


^2 

T 


> c 


dE 


(3.5) 


Several  target  fluctuation  models  were  studied  by  Swerling  [7]. 
In  particular,  we  consider  the  Swerling  II  and  IV  models. 

3.1  Swerling  IV  Targets 

If  p(E)  corresponds  to  a Swerling  IV  ("one-dominant  plus 
Rayleigh")  amplitude  distribution  [7]; 


11 


P(E)  = 


E 

e o 


(3.6) 


then 


= mean  signal  energy 


o'*  > C = 

T — 


(3.7) 


Notice  that  the  above  bound  is  twicp  as  large  as  that  for 
the  known  amplitude  and  phase  case. 


3.2  Swerling  II  Targets 

If,  on  the  other  hand,  p(E)  corresponds  to  a Swerling 
II  (Rayleigh)  amplitude  distribution  17]; 


P(E)  = 


E = mean  signal  energy 
o 


(3.8) 


then,  for  this  case  in  which  low-energy  signals  are  more  likely 
than  in  the  previous  case,  the  integral  (3.5)  for  the  bound  diverges. 


We  have  therefore  shown  that  no  finite  variance  estimator  of  (non- 


random) time-of-arrival,  which  is  also  unbiased  for  all  values  of 
signal  energy,  exists  for  Rayleigh  amplitude  distributions.  For 
the  same  amplitude  distribution,  the  ordinary  Cramer-Rao  bound 
exists,  and  is  [2]: 


(3.K 


One  arbitrary  remedy  to  this  case  is  to  only  estimate  target 
range  after  a reliable  detection  of  the  signal.  Another  method, 
which  can  be  rigorously  formulated,  is  to  estimate  the  target  range 
only  at  the  vicinity  where  the  signal  is  likely  to  appear.  This 
leads  to  the  following  case. 

3.3  Swerling  II  Targets  with  A Priori  Information 

Assuming  that  the  time-of-arrival  is  random  vith  a normal 

2 

a priori  distribution  with  variance  o''  , the  standard  Cramer-Rao 


bound  is  then; 


(3.11) 


wheru  is  the  bound  of  equation  (3.10)/  i.e./  without  a priori 
information  about  t.  Our  modified  bound  is; 


where 


(3.12) 


Figure  (3.1)  shows  the  ratio  of  the  modified  bound  (3.12)  to 
the  standard  bound  (3.11)  as  a function  of  the  parameter  Ba^.  Thus, 
the  resolution  of  the  signal  improves  relative  to  the  accuracy  of 
the  a priori  information  as  the  abscissa  increases. 

So  long  as  the  a priori  variaiice  is  nonzero,  the  modified 
bound  exists  and  exceeds  the  standard  bound  except  for  low  signal- 
to- noise  ratio  and  relatively  accurate  a priori  knowledge.  We 
therefore  infer  that  for  the  Rayleigh  channel  the  modified  bound 
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MODIFIED  BOUND 
STANDARD  BOUND 


-3-1942 


(SIGNAL  RMS  BANDWIDTH)  X (A  PRIORI  STANDARD  DEVIATION) 


Fig.  3.1.  Comparison  of  modified  and  standard  bounds  on 
variance  of  time-of -arrival  estimates  for  a Rayleigh  channel 
where  time-of -arrival  is  a random  variable. 


r 
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is  usually  tighter  than  the  standard  bound  for  estimators  of  the 
random  time-of-arrival  which  are  unbiased  for  all  signal  amplitudes. 

4.  Two-Point  Targets  - Length  Estimation 

Consider  the  radar  application  in  which  two  scattering 
centers  are  spaced  closer  th?"  signal  resolution  width.  This 

situation  occurs  when  one  attempts  to  estimate  the  length  of  a target 
with  only  two  dominant  scatterers.  Let  the  echo  from  the  ith  scat- 
terer  be  denoted  as  a^  s(t-Tj^)  . The  total  received  signal  may  be 
expressed  as 

r(t)  = a^^  s(t-T^)  + a^  s(t-X2)  + w(t)  (4.1) 

where  and  a^  are  the  complex  signal  amplitudes.  The  parameter 
to  be  estimated  is  the  time  separation  of  the  targets, 

For  fixed  values  of  ] a^ | and  la2l'  ^ conventional  way  of  evaluating 
the  Cramer-Rao  bounds  leads  to  the  following  bound  [2] . 

«o  ' 

2 

where  the  = [a^j  are  energies  in  the  signals.  This  bound  is 
obtained  by  assuming  that  the  an>l  T2  a ' nonrandom  parcimeters 
and  the  arg(aj^)  and  arg(a2)  are  random,  uniformly  distxxbuted 
on  (0.2tt),  and  applies  to  joint  estimators  of  all  four  paraumeters. 


It  does  not  include  the  effect  of  interference  between  the 
signals,  hence  is  tight  only  if  the  signals  are  well  resolved. 


To  find  the  modified  bound  one  first  obtains  the  standard 
bound  for  fixed  [a^^l  and  la2i  by  evaluating  and  inverting  a two 
by  two  Fisher  information  matrix: 


^ 26 


i-  + -i- 

E,  E. 


N, 


/ 


1+2 


p(t) 

g 2 cos  6 


<4. 3) 


where  6 » arg(a^)  - argCa^) 

.2 

p(t)  * — j Ip(t)I 

dT^ 


p (t) 


a(t+T/2)  s*  (t-x/2)  dt 


= signal  autocorrelation  fimction 


Assuming  that  0 is  uniformly  distributed  in  {0,2-rr]  and  obtaining 
the  modified  bound  by  averaging  (4.3)  over  0 yield 


(a  + 

oo2  \N_ 


p(t) 


Normalizing  the  modified  bound  by  the  conventional  bounds  i.e. 
equation  (4.2),  yields  • 


C(Ej^,E2) 

‘=o‘^l-=2> 


(4.5) 


In  order  to  evaluate  the  above  ratio  numerically,  a signal  model 
must  be  chosen.  Suppose  the  signal  auto-correlation  function  is 
Gaussian-shaped 


(T)  . 


Letting  i = 6t 


= normalized  separation 


and  evaluating  (4.5)  yield 

C(E3^,E2) 


1 - (l-»,^)^  e **■ 


One  can  obtain  an  even  tighter  bound  by  treating  all  four 
unknowns,  arg(a^)  and  arg(a2)  as  par^uueters  to  be 

estimated,  evaluating  the  four  by  four  Fisher  information  matrix, 
inverting  it,  and  tl:en  averaging  over  arg(a^)  and  arg(d2)  to  obtain 
the  modified  bound. 


This  somewhat  laborious  procedure  is  carried  out  in  Appendix 
A.  The  result  is: 


C (Ej^,E2) 


where 


(Ej^+E2l'N^ 

2E.E, 

X i. 


2n 


1 r ij^(o,T,9)  <39 

2ir  J [li)  (o,T,0)  ~ (t,t,0)  cos^e] 


(4.8) 


'i^(T^,T2»6) 


ip(T^)  + 


p{Tj^)*pfx2)  sin^  e 
[1-p^  (T2)  cos^  0] 


(4.9) 


The  above  bound  normalized  by  the  conventional  bound 
(eq.  (4.2)),  averaged  over  9,  and  using  the  Gaussian-shaped 
autocorrelation  function  is  shown  in  Figure  (4.1).  Also  shown 
is  the  ratio  of  eq.  (4.7)  which  was  more  simply  obtained  by  always 
treating  arg(a^)  and  arg(a2)  as  nuisance  parameters.  Notice  that 
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when  the  separation  is  large  (£>3",  these  three  bounds  are 
identical#  while  when  the  separation  is  small,  the  modified  bounds 
are  much  tighter,  especially  when  £ is  less  then  unity.  The  modi- 
fied bound  obtained  under  the  assumption  that  arg(a^)  and  arg(a2) 
are  to  be  estimated  (eq.  (4.8))  is  only  slightly  tighter  than  that 
obtained  under  the  assumption  that  they  are  known  to  the  estimator 
(eq.  (4.7)). 

5.  Distributed  Targets 

In  this  section,  we  turn  our  attention  to  a rather  different 
class  of  radar  returns,  namely,  the  distributed  targets.  One  very 
commonly  known  distributed  radar  echo  is  the  wake  trailing  a re-entry 
vehicle.  Such  signals  also  occur  as  a result  of  disturbances  in 
communication  media.  References  [2'J  and  [8]  contain  several  examples 
of  random  signals  in  radar  and  communication  applications.  Figure 
(5.1)  is  a simple  illustration  of  random  signals  contained  in  radar 
returns . 

Althcugh  the  basic  formulation  of  the  Cramer-Rao  bound  is  the 
same  for  various  target  configurations  (Section  2) , the  explicit 
solution  of  the  Cramer-Rao  bound  for  estimates  of  parameters  of 
a random  signal  is,  however,  usually  difficult  to  obtain.  It  was 
shown  (Reference  [2])  that  the  Cramer-Rao  bound  for  a random  signal 
mi^y  be  obtained  via  the  solution  of  a rather  complicated  integral 
equation.  Except  for  a few  special  cases,  the  solution  can  only  be 
obtained  by  using  approximations.  In  this  note,  we  consider  the 
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calculation  of  the  Cramer-Rao  bound  for  a certain  class  of  random 
signals  via  a sampled  data  approach.  The  random  signal  considered 
is  assumed  to  be  a conditional  Gaussian  process  [2]  with  known  time 
varying  mean  and  covariance  functions.  Since  many  radar  data  are 
recorded  in  sampled  form,  the  solution  presented  here  seems  to  be 
useful. 

We  present  a derivation  of  our  main  result  in  the  Section  5.1. 
A numerical  example  is  given  in  the  Section  5.2.  A deterministic 
point  target  and  a Rayleigh  target  are  limiting  cases  of  the  dis- 
tributed target.  The  derivation  which  shows  this  relation  is  pre- 
sented in  the  Appendix  B. 

5.1  Derivation 

Since  the  distributed  target  is  a very  different  class  of 
target  from  those  discussed  in  the  previous  sections,  here  we 
re-state  the  Cramer-Rao  bound  with  a slight  inconsistency  in 
notation ; 

Let  "a!' denote  the  parameter  to  be  estimated.  The  Cramer-Rao 
bound  on  the  variance  of  any  unbiased  estimator  of  a can  be  expres- 
sed as 

“j'  i = ^ ‘Va  ')  } 

(5.1) 
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2 

where  o''  = variance  of  an  unbiased  estimate  of  a 

ci 

y = measurement  process. 

A distributed  target  is  characterized  by  a complex  scattering 
function  h(t,a)  which  is  a conditional  Gaussian  random  process  with 
mean 


E(h(t,a)/a)  = Mj^(t,a) 


(5.2) 


and  covariance  function 


E[(h(t,a)  - Mj^(t,a))  (h(u,a)  - Mj^(u,a))  /a] 


(5.3) 


= H(t,a)  6(t-u) 


where  denotes  the  complex  conjugate  and  6(  ) is  the  Dirac 
delta  function.  In  addition,  we  assume  that  the  measurement 
process  Y of  eq.  (5.1)  is  at  the  output  of  a matched  filter. 

Let  s(t)  denote  the  transmitted  signal,  the  returned  signal  before 
the  receiver  is 


w (t ,a) 


s ( t-x) 


h(x,a) 


dx 


(5.4) 


The  output  of  the  matched  filter  is 


y(t,a)  = J'  w(x,a)  s* (x-t)  dx  + J n(x)  s*(x-t)dx  (5.5) 


where  n(t)  is  the  receiver  noise  with  Gaussian  distribution,  zero 
mean,  and  one  sided  spectral  density  N^.  Introducting  the  signal 
autocorrelation  function 


p(t-x)  = J's(.X'-x)  s*(x~t)dx 


(5.6) 


Equation  (5.5)  can  be  rewritten  as 


y(t,a)  = J"  p(t~x)  h(x,a)dx  + J'  n(x)  s*(x-t)dx  (5.7) 


The  output  process  y(t,a)  is  a conditionally  Gaussian  random  process 


with  mean 


E(y(t,a)/a)  = My.(t,a)  = ^ p(t-x)  (x,a)dx  (5.8) 


and  covariance  function 


S(tj^,t2,;a)  = y*  p(tj^-x)  p(t2~x)  H(x,a)dx  + N^p(t^-t2)  (5.9) 


In  many  radar  and  communication  applications,  the  data  at  the 
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matched  filter  output  are  often  recorded  and  processed  in  sampled 
form.  Let  and  denote  the  data  vectors  representing  the 

sampled  outputs  from  inphase  and  quadrature  channels  from  the  ith 
set  of  measurements  (or,  a set  of  sampled  data  of  the  ith  radar 
pulse).  Furthermore,  it  is  assumed  that  the  significant  portion 
of  the  signal  energy  is  included  in  the  data  vector.  The  covariance 

of  Y . . is 
ID 


R = 


S (t^,t^ja) 


S (t2»ti;a) 


S(tj^,t2;a)  ...  S(t^,t^;a) 

S(t2,t2?a)  . . . S(t2»tj^;a) 


S(t,t2;a).,.  S(tj^,tj^;a) 


(5.10) 


; for  all  i,j 

where  M is  the  number  of  seunples  in  the  data  vector.  The  conditional 
density  for  total  N pulses  can  be  expressed  as 

(5.11) 


P , (Y/a)  = — 

(2tt1rI) 


^ {exp -1/2  £; 


i=l  j=l 


v^here  M is  the  vector  w'ith  samples  My(t^,a)  i— 1,.  .,M.  Differentia- 
ting the  natural  log  of  respect  to  a yields 
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N 


ir  Iln  (Py/a  >1  = 


5 S S{‘"ir"y 


,T  3(R~^) 
3a 


i=l  3=i 


Squaring  (5.12)  and  taking  expectation  yield 

>2  {(-Is  ) } 


(5.12) 


(5.13) 


-II)  =]•»  © .-‘(© 

«» 

where  Tr[  ] denotes  the  trace  of  the  enclosed  matrix.  To  obtain 
the  above  results,  the  following  assumptions  and  matrix  identities 
have  been  used. 

Asstimptions 

1)  The  data  from  the  inphase  and  quadrature  channels 
are  independent 

2)  For  j^k,  y. . and  Y.,  are  independent 

X 3 XK 


Matrix  Indentities 


1)  E {(Y.  .-My)  [(Yij-My)*^  <^ij-V^}  = 

2)  E Fy.  K Y.  K y.  .1  = 2Tr[KRKR]  + Tr^[KR] 

L X]  X3  IDJ 


where  K = 


3R 


3a 


0 - a zero  vector 
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3)  TriKR]  = 0 


The  final  result  for  the  Cramer-Rao  bound  is 


(5.14) 

This  equation  is  the  main  result  for  the  distributed  target  case. 

We  make  the  following  remarks: 

1)  This  formula  is  an  approximate  solution  of  the  exact 
formulation  in  Reference  [2].  It  is  a valid  approximation  if  the 
data  vector  contains  the  significant  portion  of  the  signal  energy. 
This  requirement  can  be  easily  satisfied  in  radar  applications. 

2)  Since  the  data  vectors  are  assumed  independent,  the  above 
result  applies  to  any  unbiased  estimators  which  process  N pulses 

(or  data  sets)  incoherently. 

3)  In  most  communication  and  radar  signal  processing  appli- 
cations, the  fading  channel/fluctuating  targets  are  modeled  by  using 
a scattering  function  specified  by  Mj^(t,a)  and  H(t,a).  The  Cramer- 
Rao  bound  is  computed  by  first  evaluating  equations  (5.8)  and  (5.9) 
to  obtain  M^(t,a)  and  S(t^;tj,*a).  This  step  may  require  numerical 
integration. 

4)  The  above  result  is  explicitly  shown  for  time  domain 
data.  Due  to  the  duality  property  between  time  and  frequency,  it 
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may  be  applied  to  frequency  domain  data  (e.g. , estimating  the 
Doppler  spread  of  a fading  signal)  in  a straightforward 
manner. 

5)  This  result  can  be  extended  to  the  multiple  parameter 
case  by  computing  each  terms  of  the  Fisher  information  matrix 
(Reference  [1])  with  the  same  approach  as  above. 

5.2  Excimple 

In  this  section,  we  consider  an  example  which  applies 
to  targets  with  zero  mean  and  Gaussian- shaped  covariance  fmictions, 

i.e. , 


Mj^(t,a)  = My(t,a)  = 0.  for  all  t (5.15) 


and 


H(t,a)  = E 


/fii  a. 


■ -] 


exp  i -1/2 


(t-a) 


(5.16) 


where  E is  the  signal  energy.  Notice  that  the  parameter  to  be 
estimated,  a,  is  the  center  of  H(t,a) . Assuming  that  the  radar 
signal  is  a LFM  pulse,  the  signal  autocorrelation  function  of  a 
Hamming  weighted  LFM  pulse  can  be  approximated  by 


p (t) 


e 


(1.226B) ^ 


(5.17) 
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where  P is  the  LFM  bandwidth  in  Hz.  Notice  that  with  the  eibove 
assumptions,  equation  (5.9)  can  be  evaluated  analytically.  Define 
the  following  peak  signal-to-noise  ratio  (S/N)  at  the  receiver 
output  by  using  equation  (5«S). 


r 


S/N  = J if(x)|  H(x,a)dx 


(5.18) 


Using  the  above  assumptions,  the  Cramer-Rao  bound  normalized 

to  the  inverse  of  B /~N  is  evaluated  with  respect  to  the  target  s^'»cond 

central  moment  (a.  ) normalized  to  the  inverse  of  B where  N is  the 

h 

number  of  pulses.  The  results  for  S/N  - 0,8,  and  16  dB's  are  shown 
in  Figure  (5.2).  Notice  that  the  estimate  standard  deviation  degrades 
rapidly  for  random  targets  with  large  time-spread  (cr^^)  . when  the 
time-spread  becomes  small,  the  Cramer-Rao  bound  is  asymptotic  to 
that  for  a Rayleigh  target. 

The  above  result  can  be  extended  to  parameter  estimates  of 
a RV  w'a)te.  Results  of  Figure  5.2  correspond  to  the  time-of-arrival 
estimation  of  wakes  with  Gaussian  shaped  RCS . For  wakes  with  any 
other  shapes,  one  simply  replaces  H(t,a)  by  the  underlying  model 
corresponding  to  that  shape  where  ’h"  is  the  wake  parameter  to  be 
estimated. 

6.  Summary  and  Conclusions 

Certain  appj.ications  of  the  Cramer-Rao  bound  to  problems  in 
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radar  echo  time-of-arrival  estimation  have  been  presented.  These 
applications  encompass  point  target  echos,  two-point  target  echos, 
and  distributed  target  echos. 

With  respect  to  point  targets,  application  of  the  modified 
bound  on  the  variance  of  time-of -arrival  estimates,  valid  for 
estimators  which  are  unbiased  for  every  value  of  signal  amplitude, 
shows  the  following: 

a.  the  bound  for  Swerling  IV  targets  is  twice  that 
for  nonfluctuating  targets  having  the  same  mean 
signal  energy. 

b.  No  such  estimator  has  finite  variance  for  Swerling 
II  targets  unless  there  is  a priori  information 
regarding  the  cime-of-arrival.  When  such  information 
is  available,  the  resulting  modified  bound,  equation 

(3.12),  is  usually  tighter  than  the  standard  bound 
commonly  used. 

When  applied  to  the  problem  of  estimating  the  apparent  extent  of 
a two-point  target,  the  standard  Cramer-Rao  bound,  equation  (4.2)^ is 
independent  of  the  target  extent;  whereas  the  modified  bound, 
equation  (4. 4) , does  display  the  effect  of  interference  between  the 
two  point  scatterers  and  is  uniformly  tighter  than  the  standard 
bound.  Results  shown  in  Figure  4.1  indicate  that  the  effective 
resolution  of  a signal,  for  this  purpose,  is  about  30%  to  50%  of 
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its  reciprocal  RMS  bandwidth  in  radians/sec.  This  corresponds  to 
25%  to  40%  of  the  reciprocal  instantaneous  bandwidth  for  a LFM 
pulse. 

In  applying  the  Craraer-Rao  bound  to  estimates  of  the  parameters 
of  echos  from  distributed  targets,  a sampled-data  formulation  was 
used.  The  general  result,  equation  (5 . 14), , for  the  scalar  case  was 
applied  to  the  case  wherein  the  target  is  random  having  a scattering 
function  which  is  a Gaussian-shaped  function  of  range.  Results  of 
this  analysis,  given  in  Figure  5.2,  show  the  degradation  of  the 
time-of-arrival  estimate  as  the  time-spread  of  the  target  echo 
increases. 

Although  our  discussion  has  centered  on  the  analysis  of  accuracy 
of  time-of-arrival  estimates,  the  results  are  equally  applicable 
to  problems  in  estimating  Doppler  shift,  recognizing  the  time- 
frequency  duality.  Thus,  all  of  our  analyses  can  be  made  to  apply 
to  Doppler  shift  estimation  by  simply  interchanging  frequency  and 
time  variables.  Similarly,  the  results  are  applicable  to  angle-of- 
arrival  estimation  problems  for  a linear  aperture  by  appropriately 
reinterpreting  the  time  variable  as  of f-boresight  angle  and  the 
frequency  variable  as  distance  along  the  aperture,  measured  in 
wavelengths . 


APPENDIX  A 


Derivation  of  Equation  (4.8) 

Equation  (4.^8)  was  obtained  by  evaluating  and  inverting 

the  Fisher  Infoaraiation  ioatrix  by  treating  arg(a. } and 

X X 

arg(a2)  as  parameters.  Restating  the  received  signal 

r(t)  - aj^s(t-Tj^)  + a2s(t-T2)  + w{t)  (B.l) 


The  log  of  the  conditional  density  function  is 

Inp^  * f lt(t)  - - a2s(t-T2)  1^  dt  (B.2) 


Evaluating  partial  derivatives  of  Inp^  and  taking  expectations 
yield 


0 


^ No 


-la^l  i'(o) 


laj^l  la2lp(  ')  cosO  0 


I a 1 1 a2  I ^ ( T ) sino 


|aj_|  la2lp(i)  cos.'  -Ia2l  p(o)  j a^^l  la2l  f<  t)  sinO  0 


a^\ la2|p( t)  sinO 


a^  i 1 a.,  I p ( r)  sinO 


a,!^  ' ) COSO 


a,  1 |a-  1p(t)  cosO  1 


(B.  3) 


Since  we  are  only  interested  in  evaluating  the  bound  on  v,  one 
need  not  invert  the  complete  matrix.  Instead,  we  will  apply  the 
following  matrix  identity  to  find  the  inverse  corresponding  to 
the  upper  2x2  portion  of  F . 


Matrix  Identity 


Given 


■^1 

^121 

A = 

1 

(M 

< 

— 

,-l 

= B = 

®ll 

®12 

A 

^21 

®22 

where  A. . and  B. . are  partitions  of  A and  B,  respectively,  then 
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®11  ^^11  ~ ■^12^22  ^21^ 


Applying  the  above  identity  to  matrix  F yields 


-1, 


p _ F F -^F 

11  12  22  21 


N 


■^-P(0)- 


•2 

p (t)  sin  6 


p^ (t) cos^e 


/ \ i“p  (t)cos  e / 


|aJ|a2|fp(T)t  P±L)piiT.)sin"e\ 

^ \ 1-P^(T)COS^0  / 


COS0 


/-p(o,-  f '^>7 

\ 1-p  (t)cos  0/ 


(B.4) 


4'('"ifT2»0)  = 


/ P (t,  ) ’ ("^2^  sin^0  \ 

- (p(Tj^)+  2 2 ) 

V ^ I-P^Ct.)  cos^  0 / 


(B.5) 


Equation  (B.4)  becomes 
^11  “ ^12  ^22  ^21 


a^  I i}j  (0  , T , 0) 


N 


a^l  |a2l  i|Xt,t,0)  cos0 


a^  1 1 a2  I (t  , T , 0)  cos  0 


a2l  ij^(O,T,0) 


(B.6) 
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Its  inverse  is 


la2l^<l^(O,T,0)  -la^J  la2k(T.T,0)  cos  0 ' 

1 a ^ I ^ ( 0 , T , 0 ) 
2la^|^la2l^  (^^(0,r,d)  - i};^(t,t,0)  cos  0) 

(B.7) 

2 

The  bound  on  a is  therefore 
T 

C (E^,E2,0) 


(Ej^+E2)  >i^(0,T,9)  + 2v^~E^2 

E^E2  [i|J^(o,T,0)  - V'^('r,T,0)  cos^0  ] 


where  E^=la^|^  and  E2=la2l^  were  used.  The  modified  bound  is 
obtained  by  averaging  c'(E^,E2f0)  over  0,  i.e.. 


C (Ej^,E2) 


1 

2it 


o 


C'(E3^,E2 


0)  d0 


(B.9) 
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It  however,  can  be  shown  that 


(B.l 


tp(x,T,6)  COS0 

[ijj^(O,T,0)  - iJ;^{T/T,6)  cos^e] 


de= 


because 


^I^(T,T,  6) 

[t|;^(o,T,0)  - li)^(T,T,0)  COS^0] 

is  an  even  and  positive  function  of  0.  Using  this  property, 
equation  (B.9)  becomes 

Tjj(O,T,0)  de 

[\Jj^(O,T,0)  -lp^(r,T,d)  cos^  e] 

(B.ll 


C (E 


, , - M2)  i f 

1'  2 2 J 


This  is  equation  (4.8). 


APPENDIX  B 

Limiting  Cases  of  Distributed  Targets 


In  this  appendix,  we  prove  the  convergence  of  the  distri- 
buted targets  to  two  limiting  cases,  namely,  a target  with  known 
amplitude  and  phase  and  a Rayleigh  target. 

B.l  A Target  with  Known  Amplitude  and  Phase 

Assume  that  the  target  has  a deterministic  scattering 
function  with  known  amplitude  and  a denotes  the  time-of-arrival 
of  the  returned  signal.  This  implies 


Mj^(t,a)  = /E  6 (t-a) 


(B.l) 


and 


H(t,a)  = 0 


(B.2) 


where  E is  the  signal  energy.  Using  the  above  in  (5.8)  yields 


n 


' 1^ 


My(t,a)  P(t-a) 


(B.3) 
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Substituting  the  above  in  (5.14)  yields  the  following  result  for 
noncoherent ly  integrated  N pulses. 


C (point  target) 


1 


2E 

No 


(B.4) 


wnere 


= 

3a 


When  the  number  of  samples  in  is  large,  then 


* T -1‘ 
My  R My 


6 ^ =( signal  r.m.s.  bandwidth  in 
rad/sec) ^ (B.5) 


The  above  bound  converges  to  the  Cramer-Rao  bound  for  a known 
target  (Reference  [3]  or  Section  3 above) . 

B.2  A Rayleigh  Target 

A Rayleigh  target  satisfies  the  following  properties 

Mj^(t,a)  = 0 (B.6) 

H(t,a)  = E 6(t-a)  (B.7) 

where  E is  the  target  energy.  The  Cramer-Rao  bound  then  becomes 

] 

I 

< 
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C (Rayleigh  Target) 


_1 

N 


Tr 


(B.8) 


X • T • T 

where  R = -E  (y  Y Y Y ) 


R = r + N p 

o 


T 

r = E Y Y 
P = [p(t^-tj)] 


Y 


rP  ( t^-a)  -j 

I p(t2“a)  I 


P<V^^ 


Y 


= ^ 
da 


The  inverse  of  R can  be  expressed  as  follows  by  a simple 
application  of  the  matrix  inversion  lemma. 


R~ 


1 


T -1  -1 

(Y  P Y + N^E) 


(B.9) 
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We  will  restrict  the  number  of  samples,  M,  to  be  an  odd  integer. 

When  M is  even,  the  modification  to  the  arguments  to  follow  is 
straightforward.  From  the  definitions  above,  the  (— ^)  th  rov; 

(and  column)  of  P is  ,Y.  Other  rows/columns  of  P are  just  multiple 
shifts  of  Y.  In  addition,  the  center  element  of  Y is  unity  and  that 
of  Y is  zero.  Let  q denote  a coliimn  vector  with  all  elements 
equal  to  zero  except  the  th  element  which  is  equal  to  one. 

Then 


Y = Pq 


(B.IO) 


or 


q 


P“\ 


(B.ll) 


Using  the  above  transformation  in  (B.9)  and  after  some  manipula- 
tions one  obtains 


(B.12) 


•T. 

Y Y ) - U (y  Y + 


T.l 

Y Y )J 


(B.13) 
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j 

i 

\ 

i 

■1 

A 

\ 
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where  U is  a (MxM)  matrix  with  all  the  elements  zero  except  the 
center  element  which  is  equal  to  1/(1+N^/E) . 

Squaring  (B.13)  and  after  some  manipulations  one  obtains 

Y P Y 


1 /2«E/N^)^  \ 

(Rayleigh  Target)  = -g-  ^ 7e7K„1  "t-l  ' j 


-1 


When  the  number  of  samples  is  large,  then 


• T_-l  • 

Y p y 


(B.15) 


The  above  bound  converges  to  the  Cramer-Rao  bound  for  a Rayleigh 
target  (Reference  12]  or  Section  3 above) . 


References 


1.  H.  L.  Van  Trees,  Detection,  Estimation,  and  Modulation 
Theory,  Vol.  I (Wiley,  New  York,  1968) . 

2.  H.  L.  Van  Trees,  Detection,  Estimation,  and  Modulation 
Theory,  Vol.  Ill  (Wiley,  New  York,  1971). 

3.  C.  W.  Helstrom,  Statistical  Theory  of  Signal  Detection,  2nd 
Edition  (Pergamon  Press,  New  York,  1968) . 

4.  L.  P.  Seidraan,  "Performance  Limitations  and  Error  Calculations 
for  Parameter  Estimation,"  Proc.  IEEE  No.  5 (May  1970)  . 

5.  R.  J.  McAulay  and  E.  M.  Hofstetter,  "Barankin  Bounds  on  Para- 
meter Estimation,"  IEEE  Trans.  Inf.  Theory  IT-17,  669-676 
(November  1971) . 

6.  J.  Ziv  and  M.  Zakai,  "Some  Lower  Bounds  on  Signal  Parameter 
Estimation,"  IEEE  Trans.  Inf.  Theory  IT- 15 , 386-391  (May  1969). 

7.  P.  Swerling,  "Probability  of  Detection  for  Fluctuating 
Targets,"  Trans.  IRE,  PGIT  IT-6,  269-308  (April  1960). 


' -T*fi  yrrur ' *'"^^*T*i**^ii 


SECURITY  CLA 


m 


UNCLASSIFIED 


ATION  OF  THIS  PACE  (9lu»  Data  EuUrtd) 


Hifcturi'! 

fi^xa 


REPORT  DOCUMENTATION  PAGE 


2.  GOVT  ACCESStON  NO. 

77-111  I " I 


^|E^^*ITLE  (and  Subtit 

/Ttie  Appl 


KEAD  INSTRUCTIONS 
BEFORE  COMPLETINO  FORM 


3.  RECIPIENT’S  CATALOG  NUMBER 


TYPE  OF  REPORT  I.  PERIOD  COVERED 


/ The  Application  of  the  Cramer-Rao  Bound  I 
to  EstimateB  of  Radar  Return  Time-of- Arrival  \ 
/ for  Several  Target  Configurations  • ^ 


Chaw-Bing/::hang  MM  Robert  W.^^I 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Lincoln  Laboratory,  M.I.T. 

P.O.  Box  73 
Lexiqgton,  MA  02173 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Ballistic  Missile  Defense  Prograta  Office 
Department  of  the  Army 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 


& 


! ■ 

U.  MONITORING  AGENCY  NAME  B ADDRtSi  (if  diffaniu ^ 

!!  : 

Electronic  Systems  Divisiosf'y  •) 
Hanscom  AI%  \^/  O L 

i ^ 

i , 

< 

Bedford,  MA  01731  ' ^ 

b.  PERFORMING  ORG.  REPORT  NUMBER 

Technical  Note  1977-27 


I.  ^WTRACT  OR  GRANT  NUMBER^^; 

F19628-76-C-^'2  [ 


10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  B WOMn^T  NUMBERS 


Projecm^X363y<4D215 


23  May  V77  j i 


13.  NUMBER  OF  PACES 

52 


15.  SECURITY  CLASS,  fof  this  rspott) 

Unclassified 


ISo.  DECLASSIFICATION  DOWNCkADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  I'd/ Report;  f llj\ 

Approved  for  public  release;  distribution  unlimitear”'^  l/f  ‘. 


I7i  distribution  statement  (of  tha  abstract  sntsrsd  in  Block  70,  if  differomt  from  Rspofi) 


19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number} 


Cramer-Rao  bound 
radar  echo  estimations 
target  echos 


distributed  targets 
point  targets 
two-point  targets 


Swerling  targets 
Rayleigh  target 


20^  abstract  ('Continue  on  reverse  side  if  neceeeory  and  identify  by  block  number) 

\ 

In  this  report,  certain  applications  of  the  Cramer-Rao  bound  to  problems  in  radar  echo  titne-of-arrival 
estlmatioa  are  presented.  These  applicafions  encompass  point  target  echos,  two-point  target  echos,  and 
distributed  target  echos.  With  respect  to  point  targets  and  two-point  targets,  a modified  Cramer-Rao  bound 
is  derived  and  the  results  are  substantially  tighter  ^n  those  of  the  conventional  Cramer-Rao  bound.  In 
applying  the  Cramer-Rao  bound  to  distributed  targets,  a sampled-data  formulation  is  presented  and  a 
numerical  example  corresponding  to  a wake  parameter  estimation  problem  is  given  to  illustrate  the  results. 


UU  i JAN  73  EDITION  OF  1 NOV  65  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PACE  (Whsm  Data  Enurtil 


